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The  Frictional  Resistance  and  Boundary  Layer 
of  Flat  Plates  in  Non-Newtonian  Fluids' 


By  Foul  S.  Gronville^ 

^Boundary-layer  parameters  and  frictional  resistance  formulas  for  either  laminar  or 
turbulent  Bow  are  derived  for  Bat  plates  in  power-law  non-Newtonian  Buids.  The  results 
for  laminar  Bow  are  based  on  the  known  velocity  proBles  for  pipe  Bow,  whereas  those  for 
turbulent  Bow  are  based  on  the  application  of  similarity  taws. 


Thl  possibility  of  injecting  non-Newtonian  fluids  into 
the  boundary  layers  of  bodies  to  reduce  frictional  resist¬ 
ance  raises  two  questions:  What  is  the  frictional  resist¬ 
ance  ol  bodies  in  non-Newtonian  fluids?  Or  more  funda¬ 
mentally,  what  are  the  characteristics  of  boundary  layers, 
laminar  and  turbulent,  of  non-Newtonian  fluids? 

As  far  as  can  be  ascertained  attention  has  been  di¬ 
rected  mainly  to  the  flow  of  non-Newtonian  fluids  in 
pipes  for  chemical  engineering  applications.  The  work 
of  Metzner  and  his  associates  has  been  particularly  note- 

*  This  work  was  carried  out  at  the  David  Tavlor  Model  Basin 
under  the  Bureau  of  Ships  Fundamental  Hyaromechanics  Re¬ 
search  Program. 

*  Physicist,  David  Taylor  Model  Basin,  Nav>*  Department, 
Washington,  D.  C. 

*  The  tilde  will  be  used  over  symbols  relating  to  factors  peculiar 
to  power-law  fluids. 


worthy  in  this  respect.  Dodge  and  Metzner  [I  Y  have  ap¬ 
plied  similarity  laws  to  the  turbulent  flow  of  power-law 
non-Newtonian  fluids  in  pipes. 

The  so-called  power-law  fluids  are  those  whose  charac¬ 
teristic  stress  curves  can  be  fitted  by  straight  lines  on 
log-log  plots.  In  this  respect,  Newtonian  fluids  com¬ 
prise  a  special  case  of  power-law  fluids.  The  studies  in 
this  paper  are  confined  to  flat  plates  in  zero  pressure 
gradient  as  the  simplest  of  Ixxlies  with  boundary  layers. 
By  utilizing  the  known  velocity  profiles  of  the  laminar 
flow  of  power-law  fluids  in  pipes,  local  coefficients  of 
frictional  resistance  for  flat  plates  in  laminar  flow  follow. 

The  similarity  laws  are  applied  to  the  turbulent 
boundarj^  layers  of  power-law  fluids  on  fla*  plates  to  ob- 


*  Numbera  in  br&cketa  designate  References  at  end  of  paper. 


Nomenclature^, 


.4  =*  slope  of  logarithmic  velocity 
law 

a  —  factor  in  equation  f29) 

B},  Bi  =  intercepts  of  logarithmic  veloc¬ 
ity  law;  see  equations  (22) 
and  (23/  | 

r  =  constant  of  integration  in  equa¬ 
tion  (47) 

Cf  =*  coefficient  of  frictional  resist- 
.  ance 

—  linearization  constants  in  equa- 
tion8(60)and(64) 

D  =  frictional  resistance  or  drag 
Di,  Di  s  velocity  profile  constants  in 
equations  (36)  and  (38) 
e  »  base  of  natural  logarithms 
F  *=  outer  law  function,  equation 
(19) 


/ 

/. 


friction  factor  for  pipes 
i nner  law  f unction,  equation 
(17) 


G  ^  subscript  for  quantities  at 
junction  of  inner  and  outer 
turbulent  sublayers 
H  -  shape  parameter,  H  s  i*/0 
I\,  It  -  integrals  of  outer  law  velocity 
profiles,  equations  (39)  and 
(40) 

J I  ~  transitional  sublayer  factor, 
equation  (27) 


L  = 
In 

log  =» 
n  = 

0  = 

Pu  Pi  - 

R  = 
Bd  *■ 

ft,  - 


R«  » 


r  » 
T  * 


subscript  for  quantities  at  junc¬ 
tion  of  laminar  and  transi¬ 
tional  sublayers 
natural  logarithm  to  base  e 
common  logarithm  to  base  10 
flow-behavior  index  of  power- 
law'  fluids,  equation  (2) 
subscript  for  limit  of  overlap¬ 
ping  of  inner  and  outer  laws 
constants  in  logarithmic  resist¬ 
ance  formulas,  equations 
(61)  and  (66) 
radius  of  pipe 

Reynolds  number  of  pipe 

ftD=Cf-  2R/Mn 
Reynolds  number  based  on 
length  for  power  law  fluids, 

ft,  =  Un  x/fn 

Reynolds  number  based  on  mo¬ 
mentum  thickness  for  power- 
2_1  1 

law  fluids,  Rt  =  U"  9/p^ 
radial  distance  from  center  of 
pipe 

subscript  for  quantities  at 
pinction  of  transitional  and 
inner  turbulent  sublayers 


U  =  free-stream  velocity 
U  average  velocity  in  pipe 
u  —  tangential  velocity  in  boundary 
layer 

Ur  *  shear  velocity,  Wr  =  ( r,/p)'  '* 

X  —  distance  along  boundary'  layer 
y  —  normal  distance  from  wall 
y*  nondimensional  y  for  power- 

2_i  1 

law  fluids  y*  =  ut"  A" 
ai  *  velocity  profile  constant,  equa¬ 
tion  (35) 

/Si  —  velocity  proflle  constant,  equa¬ 
tion  (37) 

3  “  boundary-layer  thickness 
3*  »  displacement  thickness 
ff  —  boundary-layer  Reynolds  num¬ 
ber,  if  ^ 

0  ~  momentum  thickness 

*  viscosity  of  power-law  fluidc; 
see  equation  (2) 

r  “  kinematic  viscosity  of  pow’er- 
law’  fluids,  y  = 
p  “  density  of  fluid 
a  local  resistance  paramete  ,  v  == 
Uluj 

T  »  shearing  stress  in  fluid 
r,  -  shearing  stress  at  w'all 


Fig.  1  Shear  curves  of  time-independent  non-Newtonian  fluids 


tain  the  usual  boundary-layer  parameters  of  momentum 
thickness,  displacement  thickness,  and  shape  parameter. 
The  overall  coefficients  of  frictional  resistance  follow  in 
terms  of  a  local-resistance  parameter.  The  subsequent 
elimination  of  this  parameter  provides  the  more  familiar 
logarithmic  resistance-type  formulas. 

An  example  is  worked  out  for  the  frictional  resistance 
of  a  flat  plate  in  an  aqueous  solution  of  a  synthetic  poly¬ 
mer  exhibiting  a  reduction  in  resistance  for  pipe  tests  in 
the  turbulent  regime.  Here  the  values  of  the  pertinent 
boundary-layer  factors  are  deduced  from  the  pipe  data. 

Non-N«wtonion  Fluids 
De/inffion 

Non-Newtonian  fluids  [2]  are  a  class  of  viscous  or  real 
fluids  for  which  the  shearing  stress  r  is  no  longer  a  fixed 
ratio  to  the  rate  of  shear;  this  fixed  ratio  is  termed  the 
coefficient  of  viscosity  for’  Newtonian  fluids.  For  two- 
dimensional  flows  the  rate  of  shear  is  du/dy  and  then 

where  u  is  the  velocity  in  the  direction  of  the  shearing 
stress  and  y  is  the  distance  normal  to  the  direction  of  the 
shearing  stress.  Representative  curves  are  shown  sche¬ 
matically  in  Fig.  1. 

As  stated  by  Metzner  [3],  Newtonian  fluids  are:  (a) 
All  gases,  and  (6)  all  liquids  or  solutions  of  low  molecular 
weight  (i.e.,  nonpolymeric)  materials. 


0  0.2  0.4  0.4  O.t  1.0 

Rtlstivt  Diilanct,  ^ 


Fig.  2  Laminar  velocity  profiles  for  power-law  fluids 


Non-Newtonian  fluids  are: 

1  Solutions  or  melts  of  polymeric  materials  of  high 
molecular  weights, 

2  Suspensions  of  solids  in  liquids,  particularly  if  the 
solid  tends  to  swell,  solvate,  or  otherwise  associate  with 
the  liquid  phase. 

“The  distinguishing  feature  of  non-Newtonian  systems 
is  seen  to  be  that  the  colloidal  rather  than  the  molecular 
properties  are  of  significance.” 

C/aaarflcoffon 

Non-Newtonian  fluids  may  be  divided  broadly  into 
three  main  categories  [3]: 

1  Fluids  with  properties  independent  of  the  time  or 
duration  of  shear. 

2  Fluids  with  properties  dependent  on  the  time  or 
duration  of  shear. 

3  Viscoelastic  fluids  which  have  some  of  the  charac¬ 
teristics  of  solids,  such  as  elastic  recovery  from  deforma¬ 
tions. 

Most  engineering  studies  to  date  have  dealt  with  the 
first  and  simplest  category,  which  will  also  be  the  one 
considered  in  this  paper. 

The  time-independent  non-Newtonian  fluids  may  be 
subdivided  into  three  categories  whose  characteristic 
shear  curves  are  shown  in  Fig.  1 : 

1  Bingham  plastics  which  require  a  finite  shearing 
stress  to  initiate  movement.  Otherwise,  the  relationship 
between  shearing  stress  and  shear  rate  is  linear  like 
Newtonian  fluids. 

2  Pseudoplastic  fluids  for  which  the  shear  curve  is 
nonlinear  and  curves  downward.  These  include  the 
majority  of  non-Newtonian  fluids. 


2 


3  Dilatant  fluids  for  which  the  shear  curve  is  also 
nonlinear  but  curves  upward. 

^ower-iow  Fiuidt 

The  power-law  fluids  are  those  characterised  by  linear 
plots  in  log-log  coordinates  of  the  curves  of  shearing  stress 
versu.s  shear  rate,  or 


These  fluids  include  pseudoplastic  fluids  n  <  1 ,  dilatant 
fluids  n  >  1 ,  and  Newtonian  fluids  n  *  1. 

Term  n  is  the  flow-behavior  index  and  4  is  usxially 
termed  the  viscosity  or  consistency  index  (often  sym¬ 
bolized  as  K).  Here  is  called  the  power-law  fluid  vis¬ 
cosity,  since  for  New’ Ionian  fluids  n  =  1,  #4  is  the  ordinary 
coefficient  of  viscosity.  However,  it  should  be  noted 
that  the  dimensions  of  fi  depend  on  n  for  non-Newtonian 
fluids. 

The  simple  analytic  statement  of  power-law  fluids 
lends  itself  readily  to  mathematical  analysis,  as  shown 
in  this  paper. 


Lomiffior  BoviMkiry  Uiyar 
VWoeffy  Profih 

In  lieu  of  attempting  to  solve  the  equations  of  motion 
for  the  laminar  flow  of  power-law  non-Newtonian  fluids 
wherein  the  difficulties  are  compounded  by  the  nonlin¬ 
earity  of  the  shearing-stress  terms,  a  simple  expedient  is 
to  assume  that  the  known  velocity  profiles  for  pipe  flow 
[3)  hold  sufficiently  close  for  the  boundary-layer  flow  on 
flat  plates.  When  the  boundary-layer  thickness  6  is 
substituted  for  the  pipe  radius,  the  velocity  profile  be¬ 
comes 


where  u  is  the  velocity  in  the  boundary  layer  parallel  to 
the  plate. 

U  is  the  free- stream  velocity  outside  the  boundary 
layer,  and  y  is  the  distance  normal  to  the  plate.  Also 

u  =  0  at  y  «  0 


u  ^  U  at  y  «  5 


Typical  velocity  profiles  are  shown  in  Fig.  2. 
For  displacement  thickness 


h*  n 
h  2n  -h  1 

and  momentum  thickness 


6  ti(n  1) 

5  ""  (2n H-  l)(3n 


(4) 


(5) 


and  shape  parameter  H  k 

n  -b  1 

The  solution  here  for  Newtonian  fluids  «  I,  is 


(6) 


H  *  2.5  (7) 

In  contrast,  the  exact  Blasius  solution  for  flat  plates 
(4  ]  which  is  sufficiently  close  is 

H  =  2.61  (8) 

The  limiting  conditions  forn  0  are  a  constant  veloc¬ 
ity  V  for 

5*  B 

y  >  0,  y  *=  0,  -  =  0,  and  H  -  2 


The  limiting  conditions  for  n  are  a  straight-line 
velocity  profile 


^  y  ?! 


1  $ 

2'  6 


1 

6' 


and  H 


3 


Shaoring  Siratt  of  Woll 

The  local  skin  friction  or  shearing  stress  at  the  wall  r« 
is  given  by 


/du\" 


(9) 


where 


/du\  ^ 
\dy/m  ~~ 


dy 


for  y  *=  0 


From  velocity  profile,  equation  (3), 

or  for  shearing-stress  coefficient, 

-  /!L±i  V  ^ 

pV^  \  n  /  h* 

where  P  ^  p/p  wUl  be  termed  the  power-law  fluid  kine¬ 
matic  viscosity. 

Substituting for  5  from  equation  (5)  yields 
(n  -b  1)* 


(10) 


(11) 


pU^  [ 


pV*  L(2n  +  l)(3n  +  2) 


T- 

J  ft/ 


(12) 


where  R#  is  the  momentum-thickness  Reynolds  number 
for  power-law  fluids, 

ft*  ^  ,!/• 

For  Newtonian  fluids  n  =  1,  equation  (12)  becomes 

ft# 

In  contrast,  the  exact  Blasius  solution  [4]  is 

Pt/*"  ft# 


(13) 


(14) 


3 


Fig.  3  Typicftl  inner  lew  velocity  profiles  (besed  on  deu  of  reference  [!)) 


Consequently  e<iuation  (12)  will  be  arbitrarily  altered 
to  agree  with  the  exact  solution  in  the  Newtonian  ca.se,  or 


^.=  0.824  [ 


(n  -I-  1)» 

(2n  +  l)(3n  +  2). 


Turbulont  Boundory  Layor 
inner  low  or  low  of  the  Waft 

Similarity  laws  for  the  turbulent  l>oundary  layer  of 
nomNewtonian  fluids  on  flat  plates  may  be  deduced  in 
the  same  manner  as  those  for  pipe  flow  [1  ]. 

Close  to  the  wall  the  mean  velocity  u  of  the  turbulent 
flow  of  a  particular  non-Newtonian  liquid^  parallel  to  the 
wall  is  considered  to  depend  on  the  normal  distance  y 
away  from  the  W'alb  the  shearing  stress  at  the  wall, 
and  the  density  p,  the  power-law  viscosity  p,  and  the 
floAv-behavior  index  n  of  the  ftiiid,  or 

«  =  /(»,  r.,  p,  p,  n)  (16) 

By  dimensional  analysis  the  variables  can  be  grouped 
significantly  in  the  following  nondimensional  ratios: 

-  =  Mr,  n)  (17) 

Wr 


is  the  power-law  fluid  Reynolds  numlx'i*  for  the  inner 
law,  and  f  =  p[/p  is  the  kinematic  viscosity  for  power-law 
fluids.  Fig.  3  shows  characteristic  velocity  pi-ofiles  for 
power-law  fluids.  The  inner  law,  equation  (17),  reduces 
to  the  well-known  Newtonian  case  for  n  =  1. 

Outer  Law  or  VotoeHy^OaheL  Law 

At  some  distance  away  from  the  wall  and  for  the  re¬ 
mainder  of  the  bounday-layer  thickness,  the  velocity 
defect  f/-u,  for  Newtonian  fluids,  has  been  found  experi¬ 
mentally  to  be  independent  of  viscosity  p  and  is  only  a 
function  of  r«,  p,  and  distance  S-i/.  Hence  for  power-law 
fluids  the  velocity  defect  may  be  assumed  to  be  inde¬ 
pendent  of  A  and  n.  Dodge  and  Metzner  (I  ]  make  the 
velocity  defect  independent  of  fi  but  dependent  on  n. 
This  is  inconsistent  with  the  concept  of  the  outer  law  as 
pointed  out  to  the  author  by  Tulin  and  Owen  Phillips  of 
Hydronautics,  Inc. 

Then 

(/  -  u  =  /(r.,  p,  y,  S)  (18) 

or  by  dimensional  analysis. 


where  m,  =  (r»/p)*/*  is  friction  or  shear  velocity, 

*  Since  the  turbulence  level  may  be  affected  by  other  rheological 
effects  such  ae  the  development  of  normal  strewoe,  it  ia  preferable 
to  specify  the  similarity  laws  only  for  each  particular  liquid. 


The  characteristic  velocity  profile  is  shown  in  Fig.  4. 
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Fig.  4  Outer  law  velocity  profile 


LogarHhmic  Ve/ocrfy  Low 

Within  the  region  of  the  boundary  layer  where  the 
inner  and  outer  laws  overlap,  a  logarithmic  expression 
results  as  a  consequence  of  the  analytic  requirements. 
This  region  will  also  be  termed  the  inner  turbulent  sub¬ 
layer. 

Eq\iating  the  derivative  of  velocity  u  with  distance  y 
of  the  inner  and  outer  laws,  equations  (17)  and  (19), 
yields 


or 


du  _  d/i  _  _  ^  dF 

dy  dy*  d(y/S) 


dF 
6/d(y/6) 


=  A 


(21) 


Since  the  left-hand  side  of  equation  (21)  is  only  a  func¬ 
tion  of  y*  and  n  and  the  right-hand  side  of  (21)  is  only  a 
function  of  y/6  they  may  be  equated  to  a  factor  A  which 
is  independent  of  y/6,  and  n.  From  the  left-hand 
side  of  equation  (21)  there  results  after  integration, 


^  In  5**  +  B,(n)  (22) 

Wr 


and  from  the  right-hand  side, 

F  =  -A\nl+  B,  (23) 

Ur  O 


Factor  Bi  is  necessarily  a  function  of  n  from  the  integra¬ 
tion  of  a  partial  derivative.  Figs.  3  and  4  show  the 
logarithmic  velocity  profiles. 

Since  both  the  inner  and  outer  laws  hold  in  the  over¬ 
lapping  i^ion,  adding  equations  (22)  and  (23)  results  in 


—  =  i4  1ni;-|-  Bi(n)  -f-  5*  (24) 

Mr 

where 

^  pi/" 

Also 

fj  =  exp  (<r  “  Bi  -  B|)^  (25) 

Suhiayrt 

Various  sublayers  may  be  distinguished  in  the  bound¬ 
ary  layer  according  to  the  behavior  of  the  velocity  profiles, 
as  indicated  schematically  in  Fig.  5.  These  are: 

1  The  laminar  sublayer  next  to  the  wall  wherein  the 
turbulent  fluctuations  are  effectively  dampened  out. 

2  The  transitional  sublayer  wherein  the  shearing 
stresses  are  affected  by  both  laminar  and  turbulent  con¬ 
tributions. 

3  The  inner  turbulent  sublayer  wherein  the  inner  and 
outer  laws  overlap. 

4  The  outer  turbulent  sublayer  where  only  the  outer 
law  prevails. 

Velocity  Law  for  tho  tominar  Sublayor 

The  inner  law,  equation  (17),  holds  here  however, 
with  no  specification  as  to  the  exact  functional  relation¬ 
ship.  Within  the  laminar  sub-layer  for  power-law  fluids 

'-'(0 

For  the  thin  laminar  sublayer  r  —  r„  and  with  boundary 
condition  u  =  0  at  y  ^  0,  there  results 
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I  WoM 

'  Lomtnor  Subloytr  ^ 

Fig.  5  Sublayers  within  boundary  layer 
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•j|i  -  (//  6)c,h//  a),.j 
Then  HumtiuK  (2:1).  aiul  (21))  Vi’<m liters 

'  -2I„  ' 

<//  (//  6)(.  .1 

whore  (//  5),,  is  to  Ih»  (lotormiiiod  implieilly. 

I’ig.  4  sliows  a  plot  of  (lie  paraholie  v(*Ioei(y  profile. 

Boundory-loyer  Parameters 

I'he  lMmiulary-lay<*r  paiameleis  6*,  aiul  U  ^  6*  $ 
are  ohtatiu^l  from  (!»'  similarity  laws  l)y  appropriately 
integration  n  with  r<*s|XM*t  to  u  pieeiMvise  across  each 
siihlayer.  'I’he  results  are  show  ii  Table  I  w  hich  is  similar 
to  that  in  reference 
'fhen 


^ '  (/>.  - 

d  a  \  ri  / 


(:12) 


which  also  agrees  with  the  inner  law. 

Velocity  low  for  the  TransHionof  Sublayw 

The  transitional  sublayer  is  lM)iin(le<i  l)y  the  laminar 
sublayer  and  the  inner  turbulent  sublay(»r.  I' or  the 
Newtonian  ease  [5],  the  velocity  profile  for  the  transi¬ 
tional  sublayer  as  originally  derived  bv  Sijuire  is 

-  =  .4  In  {y*  -  ./,)  +  «,  rJ7) 

Wr 

where 

./,  =  B,  -  ^  In  ^  =  ft*  -  .1  (28) 

..4 

This  relation  starts  at  the  outer  edge  of  the  laminar  sul)- 
layer  and  merges  asymptotically  with  the  logarithmic 
velocity  law*,  eciuation  (22). 

It  w  ill  now  be  assumed  that  a  similar  relationship  holds 
also  for  power-law  fluids,  with  /b.  Ji,  and  yt*  iK'ing 
fuiK'tions  of  ».  Fig.  3  shows  typical  plots  of  the  velocity 
profiles  for  the  transitional  sublay (‘r. 

Vefocrfy  law  for  the  Outer  Turbulent  Sublayer 

For  Newtonian  Huids  Flama  1(>|  fitted  a  parabola  to  the 
nonlogarithmic  part  of  the  outer  law,  or 

Sr' 2)' 

A  similar  relation  will  be  assumed  for  non-Newtonian 
fluids. 

This  relation  is  to  merge  smoothly  with  the  logarithmic 
relation  equation  (2:i);  that  is,  with  e(iual  tangents. 
Hence,  if  the  point  of  junction  is  (//  6)c,  eipiating  deriva¬ 
tives  results  in 


and 

r.-- 

-K'-t)] 

(••<4) 

where 

1 

11 

a 

(35) 

~  G)<; 

(36) 

.K.vt,*- 

~  A5>i.*  +  A*) 

(37) 

+  AM  +  h 

(38) 

m 

1 

1  I 

(39) 

and 

[■  -  (f).]' 

(40) 

Since  y  is  given  in  terms  of  <r  in  ecpiation  (25),  6*,  6,  and 
H  are  then  functions  of  a  (and  n)  according  to  the  pre¬ 
ceding  relationships.  * 

Frictionpl  RttUtanca  of  Flat  Plotot 

Momenfum  fquoffon 

Considerations  of  the  momentum  changes  of  any 
fluid  [4],  Newtonian  or  non-Newtonian,  flowing  past  a 
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r 


flat  plato  iiulicute  that  the  von  Kurmun  moiiKMitum 
(Hiuutinn  is  also  applicahle  to  non- Newtonian  fluids,  or 
for  flat  plate  in  zero  pressure  gradient 


(Lv  “ 


(41) 

Then 


_  ft  ("i  _  4  ^  li 

<r  L  n  <r 

_J.  - ^(2/«)-l  ^4g)« 

2  2  —  n 


w  here  .c  is  the  stream  wise  distance  from  the  leading  edge. 

Sin<*e  the  frictional  resistance  or  drag  I)  for  a  flat  plate 
of  unit  hreadth  is 


"■r 


+  2.4  ( 

1 

tM' 

(42) 

V 

n  Di)  <r*  J 

(lie  drag  coolficiont 


+  (1  -  n)a,  + 


The  solution  of  C/  as  a  function  of  ft,  is  hence  given  im¬ 
plicitly  in  terms  of  <r  by  equations  (34)  and  (49). 

Tor  Newtonian  fluids  n  =  1,  equation  (49)  reverts  to 
(4‘^)  that  given  in  reference  (7 ). 


\\  here 


and 


iamfnar  flow 


pifH 


ft.  - 


Substituting  the  expression  for  shearing-streas  coef¬ 
ficient  for  laminar  flow,  ecpiation  (15),  into  the  momen¬ 
tum  ecpiation  (41)  and  integrating  from  .r  =  0  produces 


Byaluation  of  Constant  of  Integnthn 

The  evaluation  of  the  constant  of  integration  C  de¬ 
pends  on  the  starting  point  of  the  turbulent  flow.  If 
the  turbulent  flow  is  assumed  to  start  at  the  leading  edge, 
it  is  found  that  the  similarity  laws  do  not  hold  all  the 
way  but  begin  only  at  the  point  where  the  overlapping 
of  the  inner  and  outer  laws  starts.  The  point  is  given 
hy  yr  —  yo^  where  yr  is  the  limit  of  the  overlapping 
region  nearest  the  wall  and  yo  the  farthest.  Then 


and 


(.50) 


=  0.824*/ 


^n+l  _ 

S- 


4-  J)2  +  (l/«)  -|n/(it  +  l) 


2n  -h  I)(3n  -|-  2) 


^i/oi+i)  (44) 


The  drag  coeffi(‘ient  for  flat  plates  in  the  laminar  flow  of 
non-Newtonian  fluids  from  equation  (43)  becomes 


Cf  =  2(0.824) 


,i/(rt  +  i) 


(n  +  i)2+(i/»)  “|«/{»+i) 
,(27m)(3n  H^2)J 


1 


^^n/<n+l) 

(45) 


Turbulent  flow 

Since  s  r,c/pU^,  the  momentum  equation  (47) 
becomes 


R.  =  J  <r^dR,  (46) 

or  integrating  by  parts 

R.  =  <r%-2jR,ad<r +  C  (47) 

From  the  statement  for  R«,  equation  (34),  and  n,  equa¬ 
tion  (25)  in  terms  of  <r,  there  results  after  iirtegration  by 
parts 

J Id.  j^l  -  +  .  .  .  j  - 


<ro  =  i4  In  ^0  4“  4“  ft  (51) 

where  the  subscript  0  refers  to  the  limit  of  overlap. 

The  evaluation  of  the  constant  of  integration  for  the 
case  of  complete  turbulent  flow  then  requires  the  assump¬ 
tion  of  the  variation  of  <r*  with  R«  from  the  leading  edge 
to  the  limit  of  overlap.  As  shown  in  referen(‘e  (8]  for 
Newtonian  fluids,  the  curve  of  a-  versus  R«  starts  at  zero 
where  it  is  tangent  to  the  laminar  line  and  gradually 
reaches  the  limit  of  overlap  where  it  becomes  tangent  to 
the  curve  specified  by  the  similarity  laws,  equation  (34). 
A  like  proc^ure  may  be  applied  to  non- Newtonian  fluids. 

The  more  realistic  procedure  is  to  have  the  initial  por¬ 
tion  of  the  flat  plate  laminar  and,  hence,  with  a  constant 
of  integration  specified  for  each  particxdar  point  of 
transition.  The  constant  of  integration  is  evaluated 
from  the  momentum  thickness  developed  by  the  laminar 
flow,  or 

ft*Uin  ~  ft^tMrb  (52) 

at  the  point  of  transition.  The  value  of  <r  at  the  point 
of  transition  is  determined  implicitly  from  e(iuation  (34). 
There  the  value  of  C,  the  constant  of  integration,  is  ob- 


•  The  terms  containing  $i  in  equations  (48)  and  (49)  do  not  hold 
for  n  -  2.  Instead  for  n  *  2,  the  term  is  — /Ji  In  »  for  equation 
(48)  and  (2  In  ^  —  1 )  for  equation  (49). 
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Table  1  Velocity  Profile  Integrals  for  Pipe  Plow 


RigiM 

Lieits 

Velocity  Law 

f— (Ty* 

J  Ur 

J(; 

-p-, 

Laiiiar 

Sablayar 

_ 1 

—  -y 

«T 

I  -.2 
jyi 

1  ^.3 

JVt 

Triatitioaal 

Siblaytf 

r*T 

laaar 

Turbulent 

Sublayer 

^•'(4 

U 

—  -  it  In  y  +  C, 

«r  ‘ 

Outer 

Turbulent 

Sublayer 

— -ff-A’ 

Wt 

""{HE 

|<r*-2<r/,  +  /j|. 

Whole 

Boundary 

Layer 

f(Oj-  JirOjy-  .^1 

m'* 

■)  '’’Ln 

\o  j 

D,  .  [(Fc+^)»  +  ^>1 W, 

8 


c 


taiiiod  from  this  valiip  of  a  and  the  ft,  conpsponding  to 
of  the  laminar  flow. 


If  the  term  involving  ia  linearized  with  respect 

to  1/r//*  or 


logor/ttimfc  Aeiistonc*  formula* 

The  elimination  of  parameter  <r  from  equation  (49) 
produces  the  more  familiar  logarithmic  resistance  formu¬ 
las  with  f  V  a  function  of  R,.  The  procedure  starts  by 
combining  e<i nations  (4.‘I)  and  (47)  to  give,  after  neglect¬ 
ing  the  constant  of  integration, 

1  ('  i  ^  f 

Thp  expressions  for  J*  ftitrdff  and  ft,  from  e((uationK 

(48)  and  (:t4)  are  inserted  into  equation  (53)  to  give  with 
ai  and  01  dropped  as  negligible 


a-  2  ]_ 


2A  ^  2.4  ^(2  - 
a  n 


n)  1 


+ 


(54) 


Through  reiteration  a  is  replaced  by  the  brack- 

ets^so  that 


and  by  the  binomial  expansion 


(55) 


m) 


and  by  inversion 

'=(r[--<r 


(57) 


Now,  after  substituting  for  ^  from  equation  (25)  and 
ignoring  ai,  0i,  and  the  constant  of  integration,  Rx  in 
ecpiation  (49)  is  written  in  logarithmic  form  as 


In  Rx  = 


<T 

A 


A 


5.+  ln/>. 


1 

n 


+  ln[l-(2S  +  g)i+^  m 


Substituting  the  appropriate  expressions  for  1/a 

and  a  from  equations  (55),  (56)  and  (57),  and  expanding 
the  logarithm  as  a  series  results  in  common  logarithms  as 


IV'* 

_  _ n  _  (A+ 

2.302«v/2  \2  A/ 


log  ftr'C/ 


\/'2n  / 
2.3026.4  \ 


(A)>'>  + 


2.3026 


n(gi  +  B,) 
2.3026i4 


+  log  2£)i* 


(59) 


wherein  terms  of  higher  order  than  (C/)*'*  have  been 
neglected. 


C/'«  =  e.+^, 
then  e< (nation  [65]  l)eeome8 

log  Rx"C/  =  -f  Qj 


where 

J\  =  ^ 

y?- 

Dt\  c,  - 

2.3026 

Di)  V2.. 

and 

fl 

(Bi  +  B,)l 

^  2.3026  L  \2 

A  J 

4-  log  2/>i" 

m 


(61) 


(62) 


(63) 


I'or  n  =  ],  Newtonian  fluids,  equation  (61)  reduces  to 
the  well-known  KArmAn-Schoenherr  formula. 

Kurthermoro,  if  log  is  linearized  with  respect  to 
1  or 


log  A  =  c,  + 

(64) 

then  e((uation  (61)  becomes 

log  ftx"  =  ^  +  Q. 

(6.5) 

or 

r  - 

(66) 

1 

(log  ft,"  -  Q,)= 

where 

=  l\  -  r. 

(67) 

and 

Q,-Q>-  c. 

(68) 

Equation  (66)  provides  an  explicit  relation  between  Cf 
and  Rx. 


PipM  Flow 

Since  some  of  the  frictional  resistance  properties  of 
non-Newtonian  liquids  may  be  deduced  from  pipe  flow', 
it  is  appropriate  to  show  in  detail  the  application  of  the 
similarity  laws  to  pipe  flow.  The  friction  factor  in  pipe 
flow  is  stated  in  terms  of  average  velocity,  which  pro¬ 
vides  a  measure  of  flow-carrying  capacity  of  the  pipe. 

Average  Ve/ocfty 

The  average  velocity  of  flow  in  a  pipe  0  is  gi\'en  by  the 
rate  of  discharge  divided  by  the  cross-sectional  area  or 

(69) 
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TabI*  2  Velocity  ProMo  brtogralt 
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TabI*  2  Velocity  Profll*  Intograb 
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Fig.  6.  Comparison  of  coefficient  of  frictional  resistance  for  fiat  plates  in  turbulent  fiow 


Then  from  equations  (71)  and  (24)  there  results  for  high 
values  of  rj 


(r-'” 


V  "h  Bi{n)  -|-  2Di  +  7)3  (74) 


or  from  the  definition  of  rj  another  form  is 


2.3026^ 


log  [(/>/*)(2/")-l 


2.3026.4 

\/2 


(i+i)log2  +  ^)  +  ft-^|±®-  (75) 


~  [/(a/n)-!  2R/P^f^ 


Numerical  ExampI* 

To  illustrate  the  application  of  the  formulas  derived 
from  the  similarity  laws,  the  frictional  resistance  of  a 
fiat  plate  in  a  non-Newtonian  liquid  is  calculated  from 
results  in  pipes. 

The  example  choseti  is  a  liquid  used  in  fracturing  oper¬ 
ations  for  oil  wells  which  shows  a  lower  friction  factor 
than  untreated  water  [9] :  0. 18  percent  aqueous  solution 
of  a  synthetic  polymer  with  n  ~  0.66  and  pi  —  2.0  X 
10“*  Ib-secVsq  ft.  A  fit  of  the  pipe  data  gives 

=  4  log  [(/•'*)*  ’'  ’  +  17.38  (70) 

Comparison  with  equation  (7.'j)  then  yield.s 
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B,  +  B,  -  2D,  +  D,  =  28.00  (77) 

From  Newtonian  flow  [4| 

yii  =  4  log  (/>'««»)  -  0.40  (78) 

and  Bi  =  5.1, 

Then 

Bi  -  2Z)i  -h  D,  =  -3.11  (79) 

Then  for  the  polymeric  solution  Bi  =  31.1  from  equation 
(77).  For  both  pipe  and  flat-plate  flow  [I,  10) 

2.3026i4  =  5.(>6  (80) 

while  Bi  =  4.9,  reference  (10],  for  flat  plates  and  Bi  = 
5.1,  reference  |1],  for  pipes. 

Hence  for  the  polymeric  solution  Bi  is  taken  as  30.9 
for  flat  plates. 

From  values  of  Bi  and  Da  obtained  from  reference  [7  ] 

log  R/C,  =  -  3.108  (81) 

for  flat  plates  in  the  polymeric  solutions. 

A  comparison  of  the  drag  coefficients  for  a  flat  plate 
in  water  and  in  the  polymeric  solution  is  shown  in  Fig.  6 
for  lengths  of  10  and  100  ft,  respectively.  A  decided 
drag  reduction  is  indicated. 
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